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RESEARCH  ON  SHIP  HYDRODYNAMICS:  1982-1984 


INTRODUCTION 

In  the  interval  from  January  1982  through  December  1984,  research  in  two 
areas  of  ship  hydrodynamics  was  performed  at  the  Institute  of  Hydraulic 
Research  of  the  University  of  Iowa.  One  was  an  extended  experimental  study, 
conducted  in  the  IIHR  towing  tank,  of  the  boundary  layer  and  wake  of  a  ship 
model  at  various  Froude  numbers.  The  other  consisted  of  the  development  of 
procedures  for  calculating  more  accurately  the  potential  flow  and  the  wave 
resistance  of  sharp-edged  bodies  such  as  the  Wigley  ship  form. 

The  staff  working  on  these  problems  included  a  Postdoctoral  Visitor,  L.K. 
Forbes,  three  graduate  students,  A.  Shahshahan,  P.-P.  Hsu  and  S.  Ju,  and 
Professor  Emeritus  L.  Landweber,  who  served  as  the  Principal  Investigator. 

Since  the  detailed  content  of  our  accomplishments  is  presented  in  the 
eleven  publications  listed  at  the  end  of  this  report,  rather  than  repeat 
these,  it  will  be  more  useful  and  interesting  to  comment  on  the  reasons  for 
selecting  the  course  of  the  research,  the  significance  of  the  results,  and  to 
indicate  directions  of  future  research  on  these  problems  that  we  believe  would 
be  fruitful. 

Resistance  of  the  Wigley  Ship  Form 

In  the  first  Workshop  on  Ship  Wave  Resistance  Computations  at  DWTNSRDC  in 
1979,  it  was  found  that  computed  values  of  the  wave  resistance  of  the  Wigley 
form  could  not  be  directly  compared  with  experimental  results,  because  the 
latter  were  obtained  with  the  model  free  to  sink  and  trim;  but,  for  the  for¬ 
mer,  zero  trim  and  sinkage  were  assumed.  Consequently,  a  study  of  the  resis¬ 
tance  of  the  Wigley  form,  with  and  without  trim  and  sinkage  restraints  was 
undertaken.  It  was  found  that  the  effects  of  the  restraints  on  the  wave 
resistance  were  large  and  that  the  corrections  for  these  restraints,  intro¬ 
duced  by  the  Workshop  Committee  for  the  Wigley  form,  and  estimated  separately 
by  C.Y.  Chen  and  F.  Noblesse  (JSR  V.  27,  No  4,  1983),  were  inadequate.  This 
work  served  as  the  basis  of  the  M.S.  thesis  of  S.  Ju  and  was  published  as  an 
I IHR  report  [ 1 ] . 

*  Numbers  in  [  ]  refer  to  the  list  of  publications  at  the  end  of  this  report. 
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Boundary-Layer  Measurements  on  the  Wigley  Form 


Pressure  taps  on  the  hull,  a  five-hole  pitot  tube,  five  pressure  trans¬ 
ducers,  and  a  digital-analog  computer  were  used  to  survey  the  boundary  layer 
and  wake  at  four  Froude  numbers.  Additional  data  at  "zero"  Froude  number  were 
taken  in  the  IIHR  wind  tunnel  on  a  double  model  of  the  Wigley  form.  This 
painstaking  towing-tank  investigation  was  performed  by  A.  Shahshahan  as  part 
of  his  Ph.D.  thesis  12],  an  abbreviated  version  of  which  has  appeared  as  IIHR 
Report  No .  302  [  3 1  • 

The  original  motivation  for  this  work  was  to  obtain  measured  values  of 
the  displacement  thickness,  which  are  required  in  a  method  of  correcting  wave- 
resistance  calculations  for  viscous  effects  (4).  Shahshahan  did  apply  his 
data  for  this  prupose  and  found  that  the  correction  greatly  improved  the 
agreement  between  the  calculated  and  "measured"  wave  resistance  [2,31,  showing 
that  it  was  important  to  take  into  account  the  effect  of  the  boundary  layer 
and  wake  on  the  wavemaking.  Of  equal  interest,  however,  are  the  effects  of 
the  ship  waves  on  the  boundary  layer.  Shahshahan  discusses  this  in  his  thesis 
[2],  and  a  comparison  of  his  data  with  the  results  of  boundary-layer  calcula¬ 
tions  is  given  in  two  papers  [5,61.  The  huge  mass  of  data  recorded  by 
Shahshahan,  which  is  stored  on  tape,  can  be  furnished  to  researchers  working 
on  methods  for  computing  turbulent,  three-dimensional  boundary  layers,  for 
comparison  with  their  predictions. 

Irrotational  Flows 

1.  Ship  rolling 

The  problem  of  determining  the  added  moment-of-inertia  of  a  rolling  ship 
section  was  the  subject  of  P.-P.  Hsu's  MS  thesis  [7].  The  low-frequency  limit 
was  ejqpressed  in  terms  of  the  coefficients  of  the  Laurent  series  which  maps 
the  ship  section  into  a  circle.  Results  were  given  for  Lewis  forms  and  a 
square  section.  This  work  supplemented  a  previous  one  by  L.  Landweber  (JSR, 
Sept.  1979)  in  which  an  expression  for  the  high-frequency  limit  of  the  added 
moment-of-inertia  was  derived. 
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2.  Viscid-Inviscid  interaction 


In  this  work,  the  displacement  effect  of  a  boundary  layer  and  wake  is 
represented  by  a  source  distribution,  either  on  the  hull  surface  or  on  its 
centerplane.  It  was  presented  by  Landweber  in  the  1978  David  Taylor  Lectures 
and  in  the  1981  Weinblum  Memorial  Lectures.  The  latter  were  published  under 
the  present  contract  in  1982  141. 

The  continuation  of  the  outer  flow  as  an  irrotational  flow  field  within 
the  boundary-layer  and  wake  regions  has  been  applied  by  Kang  (University  of 
Iowa,  Ph.D.  thesis,  July  1978)  to  a  very  thin  form  (Weinblum  et  al,  TMB  Report 
840,  Nov.  1952).  He  found  that  the  calculated  wave  resistance  was  in  much 
better  agreement  with  the  residuary  resistance  when  the  centerplane  distribu¬ 
tion  was  corrected  for  displacement  effect.  An  interesting  result  was  that  a 
consistent  first-order  calculation  was  in  better  agreement  than  inconsistent 
second-order  results,  and  that  only  a  consistent  second-order  calculation 
yielded  a  superior  comparison. 

The  Wigley  hull  was  selected  for  trial  of  the  aforementioned  procedure  on 
a  more  shiplike  form.  This  was  the  original  motivation  for  undertaking  to 
measure  the  boundary-layer  and  wake  characteristics  of  the  Wigley  model.  An 
important  difference  from  the  Weinblum  form  is  that  the  Wigley  form  must  be 
considered  as  slender,  rather  than  thin.  Thus  it  was  not  evident  that  the 
previously  used  procedures,  employing  centerplane  distributions  to  correct  for 
the  displacement  effect,  were  applicable.  This  led  to  the  investigation  of 
centerplane  disstributions  for  slender  bodies  described  in  the  next  section. 
As  has  been  mentioned,  Shahshahan  has  already  applied  the  displacement  thick¬ 
ness  derived  from  his  measurements  to  calculate  the  associated  centerplane 
source  distribution  and  the  corrected  wave  resistance  for  the  Wigley  model, 
with  promising  results.  This  is  a  tentative  result,  however,  since  only 
first-order  calculations  were  performed,  and  a  thin-ship  formula  was  used  to 
obtain  the  displacement-effect  correction  as  a  centerplane  distsribution. 

3.  Slender-body  theory 

Our  approach  to  solving  slender  body  problems  differs  from  that  of  Newman 
or  Tuck  and  Von  Kerczek.  By  interpreting  the  two-dimensional  singularities, 
which  satisfy  the  2-D  Laplace  equation  and  the  slender-body  Neumann  boundary 
condition  of  a  transverse  section,  as  elements  of  3-D  singularities,  one 
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obtains  immediately  a  velocity  potential  which  vanishes  at  infinity.  The 
usual  procedure  of  matching  inner  and  outer  solutions  then  becomes  unneces¬ 
sary. 

Our  motivation  for  developing  this  alternative  approach  was  that  we 
wished  to  apply  our  method  for  correcting  ship  wave-resistance  computations 
for  viscous  effects  to  a  form  more  shiplike  than  that  of  Weinblum,  Kendrick 
and  Todd,  but  with  its  singularities  lying  on  the  centerplane.  With  the 
latter  requirement  satisfied,  one  could  immediately  use  the  formulas  of  [4]  to 
correct  the  source  distsribution  for  the  displacement  effect  of  the  boundary 
layer,  and  the  wave-resistance  calculations  would  be  greatly  simplified. 

In  order  that  a  centerplane  distribution  be  a  slender-body  solution,  each 
transverse  section  must  have  a  2-D  axial  distsribution.  A  procedure  in  which 
the  ship  section  is  mapped  conformally  into  the  unit  circle,  and  the  solution 
of  the  exterior  slender-body  Neumann  problem  is  written  as  a  complex  Poisson 
integral  in  the  circle  plane,  was  developed  for  determining  whether  an  axial 
distribution  existed.  This  existed  if  the  complex  potential  was  singularity- 
free  within  the  region  of  the  circle  which  is  the  mapping  of  half  of  the  ship 
section.  This  procedure  was  illustrated  by  application  to  an  ogival  section 
in  [7]. 

Application  of  the  aforementioned  procedure  to  the  Wigley  ship  form 
showed  that  a  centerplane  distribution  did  not  exist,  but  that,  with  a  slight 
modification  of  its  shape,  it  would  exist.  For  this  modified  Wigley  form,  the 
axial  distribution  at  each  section  could  then  be  obtained  by  analytical  con¬ 
tinuation  into  the  section  to  the  axis.  The  analytical  details  of  this  work 
and  the  resulting  centerplane  distribution  are  given  in  [8]. 

As  is  well  known,  the  internal  singularities  of  an  ellipsoid  in  a  uniform 
stream  parallel  to  its  largest  principal  axis  lie  on  its  centerplane.  We  have 
found  an  axact  solution  of  the  slender-body  boundary  condition  for  the  ellip¬ 
soid,  a  result  that  may  be  new.  For  an  ellipsoid  having  the  same  principal 
dimensions  as  the  Wigley  ship  form,  the  agreement  between  the  exact  and  the 
exact -slender-body  solutions  were  very  good.  Substitution  of  the  latter  into 
an  iteration  formula  for  solving  the  integral  equation  for  the  ellipsoid 
yielded  even  better  agreement.  This  suggests  that,  when  a  slender-body 
centerplane  distribution  has  been  found  for  a  ship  form,  it  can  be  improved  by 
iteration  in  the  integral  equataion  for  that  form.  This  work  will  be  reported 
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in  a  nearly  completed  Ph.D.  thesis  by  P.-P.  Hsu.  The  derivation  of  the  slen- 
der-body  solution  for  the  ellipsoid  is  given  in  the  Appendix. 

4.  Calculation  of  ship  wavemaking 

Dr.  L.K.  Forbes  was  a  Post-Doctoral  Visitor  at  the  Institute  of  Hydraulic 
Research  from  January  1982  to  September  1983.  During  the  first  half  of  this 
interval  he  was  supported  by  the  present  contract  on  the  development  of  numer¬ 
ical  procedures  for  calculating  ship  wavemaking,  using  the  nonlinear  boundary 
conditions.  His  thesis  on  the  2-D  problem  of  the  wavemaking  of  a  semicircular 
obstruction  on  the  bottom  of  a  channel  [9]  had  already  treated  a  free-surface 
problem  with  its  exact,  nonlinear  boundary  condition.  This  thesis  problem  was 
extended,  during  this  period,  to  incllude  the  effects  of  surface  tension  [10]. 

On  the  ship  wavemaking  problem,  he  considered  the  numerical  treatment  of 
the  integral  equation  for  a  distribution  of  Havelock  sources  over  a  ship 
hull.  The  kernel  of  this  integral  equation,  the  normal  derivative  of  the 
Havelock  Green  function,  requires  a  very  large  amount  of  computer  time,  and  he 
developed  a  number  of  closed-form  approximations  which  are  reported  in  [11]. 
Forbes  eventually  decided  to  abandon  the  Havelock  Green  function  because  of 
its  excessive  coirputer-time  requirements  and  concentrated  instead  on  methods 
using  Rankine  sources  and  the  exact  free-surface  boundary  conditions.  His 
paper  last  year  at  the  Numerical  Ship  Hydrodynamics  Conference  at  EWTNSRDC,  on 
the  wavemaking  of  a  pressure  distribution  on  the  free  surface,  indicates  the 
direction  of  his  current  research  interests.  Several  other  papers  which  he 
published  during  his  stay  at  IIHR  are  not  listed  since  they  were  supported  by 
other  contracts. 

Possibly  more  important  than  Forbes'  specific  output  is  that  we  have 
attracted  a  talented  and  productive  researcher  to  the  field  of  ship  hydrody¬ 
namics. 
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APPENDIX 


Exact  Slender-Body  Centerplane  Distribution  for  Longitudinal  Flow  about  an 

Ellipsoid  in  a  Uniform  Stream 


The  exact  centerplane  distribution  for  longitudinal  flow  about  an  ellip¬ 
soid  in  a  uniform  stream  in  the  direction  of  its  largest  principal  axis  may  be 
expressed  in  the  form  (T.H.  Havelock,  1931;  and  T.  Miloh,  1972  and  1974). 
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Here  m(x,y)  denotes  source  strength  per  unit  area  distributed  over  the  area  of 
the  fundamental  ellipse 
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U  is  the  free-stream  velocity  in  the  positive  direction  of  the  x-axis,  and  K, 
is  the  longitudinal  added-mass  coefficient,  given  by 
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where  F  and  E  denote  the  Legendre  incomplete  elliptic  integrals  of  the  first 
and  second  kinds,  and 
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At  a  transverse  section  of  the  ellipsoid,  x  -  constant,  we  obtain  the 
ellipse 
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with  the  boundary  condition  [8] 
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where  ()>  is  the  velocity  potential,  n  is  distance  in  the  direction  of  the 
outward  normal  to  the  ellipse  in  the  plane  of  the  transverse  section  and  s 
denotes  arc  length  along  the  ellipse.  This  section  is  mapped  into  the  unit 
circle  in  the  z,  -plane  by 


where  i  =  v  -1  and  c,  =  pe 
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In  the  4 -plane  ,  the  boundary  condition  (7A)  becomes 
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since,  by  (8A),  y  =  B  cos  0  along  the  ellipse.  Also  we  have 
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Hence  (9A)  becomes 
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Since  x  is  constant  for  that  section,  the  solution  is  that  for  a  source  at  the 
origin,  given  by 


<(>  = - 2  Ux  £n  p 

a 

The  axis  of  the  ellipse,  z  =  0,  is  mapped  by  the  transformation 
an  interior  concentric  circle  of  radius 


P 


o 


,b-c  >  1/2 
b+c 


with 


y  =-|-[po(B+C)  +—  ( B— C )  ]  cos0  =  (B2-  C2  ^^cosO 
o 

On  this  circle,  we  have  by  (11A) 


9<j>  _  bc_  _x 

9p  ~  D2 

a  o 


which  yields 


,9<j>.  _  , 9<f>  )  y,9z.  -  Ufr°x  csc0 

9z  z=0~  9 p  p  /l9pJp  “  2  (B+C) 

o  o  a  p 

o 

or,  applying  (6A),  (13A)  and  (14A), 

_  Ubcx  n  x_  y2  i-l/2 

^9z  z=0  2 .  2  2.1/2  u“  2  ~  2  21 

a  (b  -c  )  a  b  -c 

Hence  there  is  a  slender-body  centerplane  distribution 


mg (x,y)  =  - 


_ Ubcx _  M  x_ 

0  2,  2  2.1/2  U"  2 

2ira  (b  -c  )  a 


,-1/2 


2  2 
b-c 


extending  over  the  area  bounded  by  the  ellipse 


2 

x 


2 

a 


+ 


2 


1 


(12A) 

(8A)  into 

(13A) 

(14A) 

(15A) 


(16A) 

(17A) 

(ISA) 


9 


which  is  slightly  larger  than  that  for  the  exact  solution  (2A).  Clearly,  the 
agreement  with  the  exact  solution  would  be  improved  by  introducing  the  added- 
mass  factor  (1+k^)  in  (18A). 


10 


r  i 


i# 


